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Abstract 



We analyze nonstationary Cs atmospheric activity concentration fluc- 
tuations measured near Chernobyl after the 1986 disaster and find three new 
results: (i) the histogram of fluctuations is well described by a log-normal dis- 
tribution, (ii) there is a pronounced spectral component with period T = 1 y, 
and (iii) the fluctuations are long-range correlated. These findings allow us to 
quantify two fundamental statistical properties of the data — the probability 
distribution and the correlation properties of the time series. We interpret 
our findings as evidence that the atmospheric radionuclide resuspension pro- 
cesses are tightly coupled to the surrounding ecosystems and to large time 
scale weather patterns. 
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I. INTRODUCTION 



Chernobyl's No. 4 reactor was completely destroyed on 26 April 1986 by explosions that 
blew the roof off the reactor building and released large amounts of radioactive material 
into the environment, particularly during the first ten days. The discharge included over 
100 radionuclides, mostly short-lived, but the radioactive isotopes of iodine and cesium were 
of special radiological relevance from a human health and environmental standpoint. The 
cesiums, with lifetimes of the order of tens of years, have long-term radiological impact. Es- 
pecially important is 137 Cs, a (3- and 7-emitter with a half-life of 30.0 years (specific activity 
87 Ci/g). Its decay energy is 1.176 MeV; usually divided by 514 KeV j3 and 662 KeV 7. 
It comprises some 3-3.5% of total fission products. (It is the primary long-term 7 emitter 
hazard from fallout, and can potentially remain a hazard even for centuries.) Radioactive 
material from the reactor plant was detectable at very low levels over practically the entire 
Northern Hemisphere. Motivated by the need for reliable inhalation dose rate estimates of 
radioactive isotopes, daily air filter samples of aerosol were collected in 1987 — 1991 using 
high- volume samplers . Filters exposed to aerosols were pressed into discs and analyzed by 
7-spectrometry to find the activity concentrations of atmospheric 137 Cs. Previous analyses 
by Garger et al. revealed an exponential decrease of the absolute values of the activity 
concentration with time. The decrease was faster than expected assuming radioactive de- 
cay only, suggesting that mechanisms such as vertical migration in soil, run-off with rain, 
or melting water and snow cover may be responsible for the reduction of the activity con- 
centration. Indeed, several experts had apparently pointed out at the I.A.E.A. Experts' 
Meeting in Vienna, August 1986, that the decay rate is faster than should be observed if 
radioactivity were the only process involved . This fast decay rate was a criticism of the 
prediction of future doses by Pavlovskij which he has since apparently accepted [|J. 

It is crucial, when analyzing such data, to know how the data were collected. (For 
example, the activity in a town near the power plant will be of a different nature from data 
collected in a nearby forest, etc.) The data we analyze were taken starting in June 1987 
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with high- volume samples in the town of Pripyat, located 4 km from the nuclear power 
plant in the northwest direction. The resuspension of particulates deposited on the land 
surface following atmospheric releases of radioactive material from Unit 4 of the Chernobyl 
Nuclear Power Station in 1986 has been studied for the past six years (see, e.g., ref. |]J). 
In that region, resuspension is a result of the action of wind blowing across the terrain and 
from mechanical disturbances such as the movement of motor vehicles and the activity of 
farm equipment. In addition, the data gathered near Pripyat reflect major disturbances 
such as the demolition of buildings and the burial of the highly radioactive "red" forest. 
These activities took place prior to 1989. Enhanced resuspension may also have resulted 
from major building and forest fires within the 30 km exclusive zone and decontamination 
works. 

The fluctuations in the activity concentrations are highly nonstationary in both space 
and time. It is well known, for example, that in the summer the radioactivity is almost 
threefold higher when passing near woodlands than when in open country. This increase is 
attributed to trees picking up radioactivity, and will result in a seasonal fluctuation with a 
period of T = 1 y. 

Previous attempts to quantify the experimentally measured probability distributions of 
the concentration fluctuations were limited to statistical characteristics important for ra- 
dioecological and radiological estimations and were inconclusive because the data did not 
fit the log-normal, exponential, or gamma distributions. Hence, statistical approaches that 
apply fractal concepts to the data analysis have been useful. Hatano and Hatano, for ex- 
ample, proposed a model with fractal fluctuation of wind speed that successfully reproduces 
data of the aerosol concentration measured near Chernobyl over a decade ||. They also are 
able to reproduce the time dependence of the resuspension factor and obtain values of the 
fitting parameters that provide important information on the emission quantity and removal 
processes of nuclides from the accident. Using their model, they predicted a power law de- 
cay of the average aerosol concentration with an exponential cutoff: (C(t)) ~ e~ i//r t~ 4 / 3 , 
and temporal correlations that decay as a power law ~ t~ 2 / 3 . Here we apply a different 
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statistical description of the data and show that after removing the exponential trend in the 
data, the residual fluctuations are well described by a log-normal distribution. Specifically, 
we find that the logarithm of the concentration is a linearly decaying function of time with 
fluctuations that are characteristic of long-range power law correlated Gaussian noise, with 
temporal correlations decaying as t~ ' 3 , i.e., much slower than the model || predicts. On 
the other hand, we were able to find no indication of the power law behavior of the average 
concentration predicted in |J. 

II. METHODS AND RESULTS 

The original activity concentration data C(t) as a function of time t are shown in 
Fig. |l](a). We apply regression to the data to obtain the exponential best fit (Fig. 0(b)) 

C (t)=Aexp(-t/r), (1) 

and empirically find A = 0.588 mBq/m 3 and r = 677 days. We "detrend" the data C(t) 
by dividing it by the exponential fit Co(t) (Fig. 0(c)). We thus obtain a dimensionless 
time series that represents relative concentration fluctuations away from the exponentially 
decaying geometric mean value. Fig. |l](d) shows the logarithms u(t) = \og l0 [C(t)/C (t)] of 
these fluctuations. (We choose a logarithmic measuring unit for the same reason that acidity 
and sound intensity are measured in pH in dB respectively.) The nonstationarity (of the 
mean) in the data, apparent to the naked eye in Fig. 0(a), appears eliminated in Fig. 0(c) 
and Fig. 0(d). 

We apply histogram techniques to find the probability distribution of u{t). Fig. 0(a) 
shows that u{t) (Fig. 0(d)) is Gaussian distributed. Therefore, the signal in Fig. 0(c) is 
log-normally distributed. Since log-normal distributions often arise from underlying multi- 
plicative processes, we study correlations in u(t). 

The power spectrum of u(t) cannot be computed using the fast Fourier transform (FFT) 
algorithms because the data are missing on up to 7% of all days. We calculate the power 
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spectrum P(f) of the sequence u(t) using a Fourier transform algorithm that can correctly 
treat the "missing" data. The "Lomb normalized power spectrum" |J, P(f), often used 
when the data are unevenly sampled or missing, is shown in Figs. |2](b) and (c) on a double 
log scale. (The Lomb spectrum gives an estimate of the harmonic content of a data set 
for frequency / that is based on a linear least-squares fitting to the function A sin(27r/£) + 
B cos(27r ft) We find prominent 1 y and 2 y spectral components in the data. This 

finding suggests that the fluctuations are due in part to environmental factors. Since the 
environment has "memory" (long-range correlated behavior), we next test for long-range 
correlations in the data. The linear behavior of the log-log plot of 

p(f) ~ r p , (2) 

indicates the presence of scale invariant behavior. The scaling exponent /3 ~ 0.7 indicates 
that the data are long-range correlated, since the autocorrelation function decays in time 
with exponent 7 = 1-/5 for < /3 < 1. Note that (3 = for uncorrelated white noise and 
(3=1 corresponds to 1/ /-type noise. We have verified our finding of long-range correlations 
using other methods and find consistent results. 

III. MODEL 

We model the long-term behavior of the activity concentration fluctuations. We first 
generate a sequence w(t) of long-range correlated Gaussian noise with /3 = 0.7 (the value of 
(3 observed experimentally). We then exponentiate this sequence to obtain a second sequence 
x(t) that is log- normally distributed as well as long-range correlated. The normalization of 
x(t) is such that the mean and variance are fixed to be identical to those of the sequence u(t) 
obtained from real data. Next, we multiply the sequence x(t) by the exponential best fit func- 
tion Co(t) found empirically, thereby obtaining the model sequence Cm(^) (Fig. 0(d))- The 
model sequence Cjvr(t) and the original data C(t) are characterized by the same exponential 
trend, log-normal distribution, and power spectrum. The predicted activity concentrations 
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never rise above 10~ 2 mBq/m -3 after 1998, even for extremely large fluctuations. The pre- 
diction is encouraging from a human health and environmental perspective, since this level 
is considered to be safe. 



IV. DISCUSSION 

The log-normal distribution || could be related || to local resuspension of the radionu- 
clides by wind and transport from distant areas. Additional support for the possibility that 
local resuspension is the leading factor determining the observed fluctuation arises (i) from 
the presence of the pronounced 1 y spectral component and (ii) from the observed long-range 
correlations, as we next discuss: 

(i) A pronounced 1-year periodicity can be attributed to large-scale seasonal environmen- 
tal changes, such as snow melting, as well as to human agricultural activity, such as 
plowing. Both yearly processes can release radioactive elements accumulated in soil. 

(ii) Long-range correlations with shorter characteristic times can be generated by large- 
scale atmospheric turbulence. It is possible that resuspension of radionuclides due to 
changing weather conditions give rise to long-range correlated fluctuations around 
the exponentially decaying average activity concentrations, Co(t). However, it is clear, 
that it is not only weather that is responsible for the observed long-range correlations 
in the activity levels. Ecosystems consist of a large number of reservoirs (e.g., forests, 
lakes, swamps and agricultural fields) of various sizes that can repeatedly absorb and 
release radionuclides into the environment [|10j. It is well known since the works of 



Hurst [11] that the level of water in lakes is characterized by power law decaying 



fluctuations with the Hurst exponent H 0.7 [11]], The amount of radionuclides 
released in a unit of time should be proportional to the derivative of these fluctuations 
and hence should have a power spectrum scaling exponent (3 = 2H — 1 [T2j , which is 



in qualitative agreement with our measurements. The value of (3 ~ 0.7 found by us is 
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larger than the value 0.4 obtained from Hurst measurements and the discrepancy may 
be partially related to the fact that the high frequency part of the power spectrum (see 
Fig. 0(c)) has smaller slope (3 ~ 0.5. Similar values of (3 ~ 0.5 as well as log- normal 
distributions of fluctuations were obtained for river discharges ||13| . Moreover, the 
radioactive release is related not only to lakes but also to other parts of the ecosystem, 
which may have different Hurst exponents. Indeed, the influence of weather (and 
possibly forest fires) on monthly mean 137 Cs concentrations provides an illustration of 
how the ecosystem can absorb and release radionuclides [ 14[] . The power law behavior 



of the power spectrum indicates that the temporal correlations decay as a power law 



~ t~ 7 , with 7 = 1-/9 JT5J. Hence our analysis yields 7 = 0.3. The model of 
wind resuspension used in || also predicted power law decaying correlations but with 
7 = 2/3, indicating a much faster decay than we observe. Indeed, such slowing down 
of the correlation decay with respect to the pure atmospheric activity is often observed 



in the the response of ecosystems such as forests and ground waters [p~3,inj] 



V. CONCLUSION 

In summary, we have attempted to find a complete statistical description of the data. 
We have succeeded in quantifying both the probability distribution as well as the correlation 
properties of C(t) — the two most fundamental statistical properties of time series. As with 
refs. fl5,|i~7l, the values obtained for the fitting parameters provide important information on 



the emission quantity and removal processes of radionuclides from the accident site. The 
differences between our results and those of ]5],[17j may indicated the significant role of the 



ecosystems in the radionuclides resuspension. One important conclusion is that the observed 
long-range correlations are not solely due to the atmospheric turbulence but are strengthened 
by the accumulation of radionuclides by ecosystems, such as forests and ground waters. 
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FIG. 1. (a) 137 Cs activity concentration C(t) measured after the Chernobyl disaster as a func- 
tion of time t in days, (b) Semi- log plot of the data and the exponential best fit Co(t) ~ exp[— 1/677], 
obtained by regression (solid line) . (c) Fluctuations after correcting for the exponential decay. The 
data C(t) have been divided by the exponential fit Co(t). (d) Logarithms u(t) = log 10 [C/Cb] of 
these data. The nonstationarity apparent to the naked eye in (a) appears eliminated in (c) and (d). 
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FIG. 2. (a) Semi-log plot of the histogram of u(t) along with a Gaussian fit, showing that u 
closely follows a Gaussian distribution (a parabola in semi-log plot). The fluctuations in C(t) are 
thus seen to follow a log-normal distribution, after being corrected for the observed exponential 
decay, (b) Normalized power spectrum P(f) of u(t) as a function of the period 1// in years. 
Due to the lack of data for up to 7% of all days, the spectrum shown has been computed using a 
Fourier transform algorithm that is immune to the "missing data" problem as discussed in the text. 
The pronounced 1 y and 2 y spectral components suggest that the fluctuations are influenced by 
environmental factors, (c) Double log plot of P(f) as a function of frequency measured in years -1 . 
We have applied logarithmic binning to smooth the spectrum by averaging over windows which 
double in size. We find that (3 (s=s 0.7) is significantly different from zero, so we conclude that the 
activity data C(t) are long-range correlated and log-normally distributed around the exponentially 
decaying function Co(t). (d) Semi-log plot of activity concentrations as a function of time after 
the nuclear accident, up to 20 years. Experimental data are shown for all the dates when they are 
available. For future dates when data are not available, we use the model discussed in the text to 
predict the long-term behavior. The predicted radiation levels for the future (12 years after the 
accident onwards) is below the generally accepted safe level. 
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